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1.  Introduction 


There  is  no  exact  test  for  testing  the  significance  of  ,  the  nesting 
effect's  variance  component  in  an  unbalanced  random  two-fold  nested  classifi¬ 
cation  model.  There  are,  however,  approximate  tests  which  utilize  ratios  of 
mean  squares  (see  Cummings  and  Caylor  1974,  Tietjen  1974,  and  Tietjen  and 
Moore  1968).  In  general,  the  exact  distributions  associated  with  these  tests 
are  complicated  which  makes  it  difficult  to  adequately  determine  the  tests' 
true  levels  of  significance  and  power  values.  A  comparison  of  four 
approximate  tests  concerning  was  recently  made  by  Tan  and  Cheng  (1984). 

In  this  paper  we  present  an  exact  F-test  for  testing  the  null  hypothesis 

=  0  versus  ^  0.  (1.1) 

Aside  from  the  usual  assumptions  concerning  the  random  effects  in  the  model, 
the  only  other  condition  for  the  validity  of  the  test  is  that  N  >  2b  -  1, 
where  N  is  the  total  number  of  observations  and  b  is  the  total  number  of 
levels  of  the  nested  factor.  The  proposed  tost  is  compared  against  the 
approximate  tests  described  in  Tan  and  Cheng  (1984)  with  respect  to  power. 

The  results  of  this  comparison  indicate  that  the  former  test  is  quite 
efficient. 

2.  The  tlevelopment  of  the  exact  test 
Consider  the  unbalanced  two-fold  nested  model 

Yj  41.  =  M  +  a.  +  8,  .  t  E.  (2.1) 

1  ij  Jjk 

1  =  1,2,. ..,a;  j  =  l,?,...,b^;  k  =  l,2,...,n^j,  where  p  Is  an  unknown  constant 
parameter,  and  8..  are  random  effects  associated  with  the  nesting  factor 
and  the  nested  factor,  respectively,  and  random  error.  We  assume 

that  and  ate  independently  distributed  as  N(0,a^),  N(0,ag),  and 

N(n,a^),  respectively.  We  also  assume  that 


iT’ 


N  >  2b  -  1, 


(2.2) 


where  N  =  "ij’  ^  ”  i-1  need  for  inequality  (2.2)  will  be  seen 

later.  We  note  that  the  latter  assumption  is  quite  reasonable  and  can,  for 

example,  be  satisfied  if  ^2  for  all  l,j. 
n  . 

Let  “  1,2,..., a;  j  =  l,2,...,bj).  From  (2.1)  we 


yjj  -  M  +  «!  + 


(2.3) 


1  =  1,2...., a;  i  =  l,2,...,b. ,  where  e . .  =  ,  c. /n. ..  Model  (2.3)  can  he 
*  ’  ’1’  ij  k=l  Ijk  ij 


rewritten  in  the  matrix  form 


V  =■  p  1.  +  A  a  +  16  +  £, 


(2.4) 


where  y  and  £  are  vectors  consisting  of  the  and  the  respectively, 

a  =  (a^a^  , . . .  ,a^)' ,  6  =  n  ♦  ^^2  ’  *  *  *  ’  ^ab  *  -^b  ®  vector  of  ones  of 

a 

dimension  b,  I,  is  the  identity  matrix  of  order  bxb,  and  A,  = 

’  -'b  ■' 

Diagfl,  ,1.  ,...,1,  ^  is  a  block-diagonal  matrix  of  vectors  of  ones.  From 

'“b,  ‘'b_  '^b 

12a  _ 

(2.4)  the  variance-covariance  matrix  of  y  is 


where 


ir  y  =  A,0^  +  I.Od  ^  Ko^, 
^  ''la  ''b  6  c' 


A,  =  A, a;  =  Q.  J,  , 

~1  1=1  'h^ 

V  n-  r  -1  -1  -1  ) 

K  =  Diagin^^,  0^2. •••  •, 


(2'-‘j) 


(2-6) 


(:''7) 


where  is  a  matrix  of  ones  of  order  b^xb^  (i  =■  1,2,..., a),  and  @  clenot 

the  direct  sun.  The  residual  sura  of  squares  for  the  model  in  (2.1)  is 

T  =  .Z.  , y, It  is  known  that  T/o^  has  the  chi-squared  distribution 
i,.l,k  ■'ijk  "'Ij  e 

with  M  -  b  degrees  of  freedom  independently  of  y.  We  can  write  T  as 

T  ^  y'  R  y,  (2.«) 


.  ■c. 


where  y  is  the  vector  of  N  observations  and  ^  is  the  matrix 


R  =  Im  -  ©•  /n.  .  V 

-  -N  fVj  ij 


(2.9) 


We  note  that  JJ  is  symmetric  and  idempotent  of  rank  N  -  b,  and  by  assumption, 
N  >  2b  -  1,  that  is,  N  -  b  >  b  -  1.  We  can  thus  express  as 


R  =  C  A  C', 


(2.10) 


where  Q  is  an  orthogonal  matrix  and  A  is  a  diagonal  matrix  whose  first  N  -  b 
diagonal  elements  are  equal  to  unity  and  the  remaining  b  elements  are  equal  to 
zero.  Furthermore,  we  can  partition  £  and  A  as 


where 


A  =  Diagfl^^,l^^,0^, 
C  =  [Ci:C,:C3], 


=  b  -  1, 

=  N  -  2b  +  1, 


(2.11) 


(2.12) 


(2.13) 


and  CpC^.C^  are  of  orders  NxUj ,  Nxv2,  and  Nxb,  respectively.  Note  that 

c:c.  =  I,  i  -  1,2,3, 

c;c.  -  0.  if) 

Formula  (2.10)  can  then  be  rewritten  as 


(2.U) 


R  =  CjC;  +  c^c',  (2.15) 

which  results  in  a  partitioning  of  the  residual  sum  of  squares  T  into  and 


To ,  wliere 


T ,  =  V '  C  C '  V  , 

1  ,V  o  Jo  I  O 


(2.16) 


t2  =  ^2^2  X- 


(2.17) 


Consider  now  the  matrix  Aj  in  (2.6).  There  exists  an  orthogonal 
matrix  P  of  order  bxb  such  that  P  A,  P^  =  A. .  where 


2.18 


where  D  is  a  zero  matrix  of  order  (b-a)x(b-a).  This  is  because  A,  has  Che 

''1 

eigenvalues  bj ,b2 , • • • .b^,  and  0  of  multiplicity  b-a.  The  first  a  rows  of 
are  orthonormal  eigenvectors  of  Aj  that  correspond  to  bj ,b2 , • • • .b^.  Let  be 


the  axb  matrix  consisting  of  these  rows,  that  is, 


P,  .  DlagQ'  //b,  .1;  //bj . IJ  //b^ 


(2.19) 


Let  ^2  (b-a)xb  matrix  consisting  of  the  remaining  b-a  rows  of  Then 

P  =  If  z  =  P  y,  then  from  (2.5)  we  have 


Var  z  =  A.o^  +  1^0^  +  P  K  P'o^, 
-la  -b3  ---e 


(2.20) 


Theorem  2.1.  There  exists  an  orthogonal  matrix  Q  of  order  bxb  such  tliat  the 
first  row  of  Q  P  Is  l'//b. 

Proof .  Define  the  unit  vector  Cj  =  (Cj:  O']  where  c'  =  (/b^j  ,/bj, , . . .  ,/b^)/i/b 

and  0'  is  a  zero  vector  of  dimension  b-a.  Then  (I,  -  e,e')e,  =  0.  The  matrix 

''b  1''  1  I  - 

I  -  e.e'  is  idempotent  of  rank  b-1.  Let  Q,  denote  a  matrix  of  order  bx(b--l) 
and  rank  b-1  whose  columns  are  obtained  via  a  Cram-Schmidt  orthonornia  ll/.at  ion 
of  the  columns  of  2  ~  0  is  an  orthogonal  matrix 

and  the  first  row  of  Q  P,  namely  e'P,  is  l'//b. 

,v  -  ''I''  ''[> 

From  Theorem  2.1  we  conclude  that  if  u  -  QTz,  where  z  P  y  and  0,  is  Che 

^  /V  /V#*.  Kl 


matrix  described  in  the  proof  of  the  theorem,  tlien 


since  Q  P  is  orthogonal,  and 


Var  u  -  Q:A,Q,o^  +  k  , P'Q.o^. 
^  -'l-'l-'l  a  -b-1  3  -1^  ^ 


(2.21) 


(2.22) 


since  QTO,  =1.  We  note  that  0'A,Q,  is  of  rank  a-1.  To  show  this  wo 

-1-1  -b-1  -l-l-^l 

partition  as  f2ii’2i2^’  2ji  (b-l)xa  and  i®  (b-l)x(b-a).  Then 

g'A.Q.  p'  Diagfb  j  ,b  , .  .  .  j  (see  2.18).  Hence,  rank  (PjAjPj)  = 


rank  ^  ^1~P  ""  ^''^’  where  c'  is  the  unit 

vector  described  in  the  proof  of  Theorem  2.1.  The  one  before  last  equality 
follows  because  the  columns  of  the  matrix  orthonormal.  It 

follows  that  there  exists  an  orthogonal  matrix  5  of  order  (b-l)x(b-l)  such 
that 

=  S  DiagfD.O)  S',  (2.23) 

where  p  is  an  (a-l)x(a-l)  diagonal  matrix  of  nonzero  eigenvalues  of 
Q  is  a  zero  matrix  of  order  (b-a)x(b-a)- 

Consider  the  vector  m  =  S'u.  From  (2.21),  (2.22),  and  (2.23)  it  can  be 
seen  that  jo  has  a  zero  mean  and  a  variance-covariance  matrix  given  by 

Var  m  =  Dlagfo.O'lo^  +  I  +  Lo^,  (2.24) 

^  a  '^b-l  B  ~  e 

where  L  =  S'Q'P  K  P^Q-S.  Define  the  vector  12  as 

s  ■  'Wh-l  -  i>'‘  Sj!!’  <2.25) 

where  is  the  largest  eigenvalue  of  the  symmetric  matrix  J.  and  is  the 

Nx(b-l)  matrix  in  (2.12).  Note  that  the  matrix  t  is  positive 

scmide  f  ini  t  e ,  hence  well  defined  with  eigenvalues  equal  to 

the  square  roots  of  the  eigenvalues  of  Let  J[2  be  partitioned 


as  n 

=  [fi':  12'] 

' ,  where 

12  and  12  „ 

-'a  -P 

-a  -P 

Then 

ren  .?  .  2  . 

(i) 

-  Q’ 

(ii) 

9^  and  12p 

are  stat 

i  s  t  i  c  a  1 1  y 

the  follciwing  va rianco-cova ri ance  matrices: 


Var  >2  =  1)0^  +  ( +  A  o^)  T 

''a  a  P  max  r  -a-1’ 


Var  n.  =  (o^  +  A  o^)  I  . 
’'P  P  max  r  ''h-a 


(2.26) 

(2.27) 


rKrwr%f§ 


(i)  This  is  true  because  K(w)  =  S  E(u)  =  0  by  (2.21),  and  =  0. 

(ii)  Since  w  =  S"  Q'  P  then  is  a  linear  function  of  the  vector  of 

observations,  hence  it  is  normally  distributed.  We  now  claim  that  and 

C'jf  in  (2.25)  are  statistically  independent.  To  show  this  we  write  ^  in  the 

form  y  =  G  y  where  C  =  Dlagfl"'  /n,  .\'  /n,„,...,K  /n  ,  ).  Since  y  is 

^  '"n,,  li’-^n,-  12  ’-n  ,  ab  * 

1 1  12  ab  a 

a 

statistically  independent  of  the  residual  sum  of  squares,  T,  in  (2.8),  then 
G  E  R  =  0,  where  E  =  Var  ^  (see  Searle  1971,  p.  59).  Using  the  representation 
(2.15)  for  g  we  obtain 


G  E  (C^C'  +  C2CP  0. 


(2.28) 


If  we  multiply  (2.28)  on  the  right  by  and  note  (2.14)  we  get 


G  E  C,  =  0. 
-^1 


(2.29) 


From  (2. "9)  it  follows  that  Cov(u  ,  J['*Cj)  =  Cov(S'gj‘  P  y  ,  y'cp  - 

rQi'  P  Cov(y.y'cp  =  S^g^P  Cov(G  y  ,  =  S-gp>  G  E  G,  0.  Hence,  the 

variance-covariance  matrix  of  g  in  (2.25)  is  of  the  form 

Var  n  Var  w  +  (X  hfi  CT  E  ('  (A  1.  ,  -  l.f-i  .  (2. 

''  max-'b-l  .vj  .v]  inax'b'l 


(2.30) 


We  claim  that  C'  E  C,  =  I, 

-1  ^1  e  ''b-1 


This  follows  from  the  f.ict  Chat  T/ of  has  the 


chi-sqviared  distribution,  hence  E  J  E  5  E  (-‘'Oe  2,B  and  Theorem  2  in 
Searle  1971,  p.  57),  whicli  can  also  be  written  as  E  E  E  '  R.  By  noting 
(2.15)  we  get 


(c.cr  +  c.,c:)  E  (c.cr  +  c.,c:)  -  (c  c:  c  g.-g,-;). 

•'I -'I  ''2-' 2  ''  -'I''!  -'2'' 2  t  "2'-2 


(2.31) 


If  we  now  multiply  (2.31)  on  the  left  by  C'  and  on  Che  vighi:  by  C,  and  note 

'^1  ''I 

(2.14)  we  obtain  the  desired  result. 

From  (2.24)  and  (2.30)  we  conclude  that 


Var  ii  =  Uiagrn,o1  +  I.  .0^  f  T,  1  (A  I.  ,  -  L)  0' 

‘  a  ''b-1  B  ''  t  inax-'b  I  £ 


(2.32) 


=  Diagfo.O] 


o'- 

a 


+  X 


max 


ib-r 


Since  Var  n  is  a  diagonal  matrix,  and  must  be  statistically  Independent. 
Furthermore,  from  (2.32)  it  can  be  concluded  that  these  random  vectors  have 
the  variance  structure  described  in  (2.26)  and  (2.27),  respectively. 

From  Theorem  2.2  we  can  then  state  that 


n'  (Do^  + 
a 


‘  la-l  '  xLl 

Ss  "e''  '  4-a' 


where  6  =  o?  A  o^.  Under  Hn:o£  =  0,  Q'Q  /<5  ~  i  ♦  hence  F  =  MS„/MSq  has 
p  max  E  u  a  ’  a-1  ot  p 

the  central  F-distribution  with  a-1  and  b-a  degrees  of  freedom,  where  MS^ 

=  n'il  /(a-1)  and  MSo  =  /(b-a).  It  is  easy  to  verify  that 

'-or-a  P''8''8  ^ 


a-1 

E(MS  1  =  6  +  [.f  d./(a-l)]o2 

a  ‘■1  =  11  ■‘a 


F,(MSg)  =  6,  (2.33) 

where  d.  is  the  diagonal  element  of  p(  i=  1 , 2  , . .  .  ,a-l  ) .  Hence,  large  value.'- 

of  the  test  statistic  F  are  significant. 

3.  A  comparison  oF  the  exact  test  against  the  approximate  tests  in  Tan  and 


i 


Cheng  (1984) 

Tan  and  Cheng  (1984)  compared  the  powers  of  four  approximate  tests  of  the 
hypothesis  Hq  described  in  (1.1).  The  corresponding  test  statistics  were 
denoted  by  Fj,  F2,  F -j ,  and  F^ ,  respectively,  and  their  power  values  were 
derived  approxim.Uely  by  using  hagiierre  polynomial  expansions  of  the  true  null 
and  non-null  di  st  ri  but  i  on.s  of  F.  (  i  1  , 2 , 3 , 4  ) .  The  power  values  were  obtained 
for  different  values  of  several  combinations  of  the  nuisance 

parameters  and  a^(-l,2,3),  and  for  two  unbalanced  nested  designs  winch  v;e 
reproduce  in  Table  1. 


1?^ 

'v 


In  Tables  2  and  3  we  give  a  listing  of  tbe  power  values  associated  with 
Fi,F2,F3,and  F^  as  reported  in  Tan  and  Cheng  (1984,  Table  4,  pp.  197-198)  for 
a  level  of  significance  a  =  .10, 


At  tbe  a-level  of  significance,  the  power  function  for  the  exact  test 
proposed  in  Section  2  is  given  by 


P(MS„/HSj  >  H,), 


where  H  is  the  alternative  hypothesis  in  (1.1),  and  F  .  , 

3  Cly3"'lji)3 


(3.1) 

denotes  the 


upper  a°<  point  of  the  F-distribution  with  a-1  and  b-a  degrees  of  freedom. 
Under  II  ,  MS^  and  MSq  are  independently  distributed  and  (b-a)MSQ/6  is 

3  CX  p  P 

distributed  as  Xk_  «  (a-l)MS  /6  no  longer  has  the  chi-squared  dlstribu- 

D  3  u 


tion.  In  this  case,  since  is  normally  distributed  with  a  zero  mean  and  a 

~a 


variance-covariance  matrix  Var  0  =  Do^  +  5l  .  ,  (a-l)MS  is 

,  -'a  a  ''a-1  ■'ct  -'a  a 

a-1 


distributed  as  >  where  the  W^'s  are  independent  chi-squared  variates 


■•th 


with  one  degree  of  freedom,  and  A.  is  the  i  eigenvalue  of  Var  ,  that  is, 

1 


th 


Aj  =  +  5  with  d^  being  the  i  diagonal  element  of  U(i -•  1 , 2  , . . .  ,a-l )  (see 


Johnson  and  Kotz  1970,  p.  ISl).  Thus,  under  H  the  exact  lest  statistic,  K  can 


be  written  as 


a-1 


F  = 


Q'  Q 

-'a 


ih  h  h 


6(a-l)MS./(5  ■  6(a-l)MS./>S 
P  P 


(3.2) 


Approximte  values  ot  the  power  function  in  (3.1)  can  be  conveniently  obtained 
by  usiiur  IMrotsu's  (1979,  pp.  578-519)  approximation  of  the  upper  probability 
values  of  a  statistic  of  the  form 

X  'A  X  /  ( c  f  ) 


H  = 


0^7  o^ 


(3.3) 


whe.c  x  is  normally  distributed  witli  a  zero  mean  and  a  variance-covariance 

,2/,,?  I,,  a ; .  t...  ►  ..,1  ci/f 


latrix  y,  A  is  a  nouiu'gative  matrix,  Is  distributed  as  (l/f2)x/. 


independeut  ly  of  x'A  x,  and  c  and  f  are  g.ivt'u  hy 


C  =  -^  x)/Kj(x"A  x) 

(3.4) 


f  =  2  <?  (x'A  x)/i<,,(x'A  x)  , 

where  k. (x'A  x)  denotes  the  cumulant  of  x'A  x.  For  convenience  we 

1  ys.  ^  ^  /«.  /s. 

reproduce  the  formula  fcr  P(H  >  h)  as  given  in  Hirotsu  (1979,  formula  2.4): 


P(H  >  h)  ■=  P(F^  f  >  h) 


[A/{3(f+2)(f+4)B(i  f,  I  f2)}](l+fh/f,) 


-  J(f^f  J 


1. 


2(f+f  )(f+4)  (fff  +2)(f+f  ) 
(fh/f  )^  [(f+2)/(f+4)  +  .11 


l+f2/(fh) 


”l+f , 


(3. 


where  ^  denotes  the  F-dist  r  ibut  ion  with  f  and  f2  dep.rees  of  freedom, 
B(mpiTi.^)  denotes  the  beta  function,  and 


k-  (x'A  x) 

t  ^  /V  />. 


<3(x 


'A  x)/k^(x'A  x) ]  -  I. 


(3. 


The  approximation  described  in  (3.5)  was  developed  via  a  I.aguerrc  polynomial 
expansion  of  the  true  distribution  of  the  statistic  x'A  x/(2c)  and  was 
reported  in  Hirotsu  (1979)  to  he  quite  satisfactory. 

From  (3.1)  and  (3.2),  the  power  function  for  the  exact  test  statistic  F 
can  be  written  as 


n'n  /(cf)  ,, 

-a-' a  ,  n(a-I  ) 


- X  p  j„  1 

MS^/5  cf  a,a-l,b-a'  a  ’ 

p 


(3. 


1  s 


tr(Var  ft  )  a-1 

''O 


tr 


;  in  (  ;i ,  4  )  ,  but 

with  it'll  subs  L  i  t.utfd 
''  ci'^  a 

a-1 

a-  I 

.i:,(d.o^  +  6)^- 
1  =  1  1  a 

5.1,(d.O  +  1)2 

1  =  1  1 

a- 1 

a  -  1 

.F,(d.a2  +  6) 

1  =  1  1  a 

.):,(d.o  +  i) 

1-1  1 

a-1 

a  -  1 

\  .lAd  ,0-  +  ^)] 
‘  1  =  1  i  a 

"  lji;,(d.o+i)}2 

a~  I 

^r-i 

(3.8) 


.l,(d.o2  +  6)^ 


1, (d. n+1 )• 


f 


(3. 


w 

>  vT 


'-'.N 

S'.N 

'•\‘- 

S'- 

I 


W: 


J.iJ  J..»..L>.L..L  ILI.U..LJ...I  u.  y  JU 


where 


0  = 


a 


=  ■;7z 


of  +  X  o'^ 
&  max  e 


(3.10) 


(see  Hirotsu  1979,  p.  579).  We  note  that  the  statistic  which  appears  in  (3.7) 

is  of  the  same  form  as  the  statistic  H  in  (3.3).  This  can  be  seen  by  taking 

X  =  ,  A  =  I  ,,  and  o^/o^  =  MS./5,  which  is  distributed  as  (l/fo)xl  with  f^ 

"'a  ''  ''a-i  p  ^  ^ 

=  b-a  degrees  of  freedom.  It  follows  that  the  power  value  in  (3.7)  can  be 
approximately  computed  by  applying  formula  (3.5)  and  remembering  that 
fl'n  /(cf) 


H  = 


MSg/5 


.  _  5(a-I )  ^  *^*^q.a-l.b-a 

cf  a,a-l , b-a 


a-1 


(using  formulas  3.8  and  3.9), 


.1,(8. 0+1) 

1  =  1  1 


■2  = 


c  and  f  are  as  given  in  (3.8)  and  (3.9),  and  A  is  as  described  in  (3.6), 

which  can  also  be  written  as 

[tr(Var  )]ftr{(Var  )^}] 

.  ‘ ^  ‘'a _ . 

[tr|(Var 


a-1  a-1 

[^|j(d.0+l)][^|j(d.0+l)a] 


[^|J(d.0+l)2]2 


1, 


(3.11) 


where  0  is  given  in  (3.10)  (see  Hirotsu  1979,  p.  579).  It  is  interesting  to 
note  that  in  (3.5)  the  power  of  the  exact  test  depends  on  and  o^ 

through  0.  A  complete  determination  of  this  power  requires  finding  the  values 
of  X  ,  d,,d«,...,d  ,,  which  depend  on  the  design  used,  and  a  specification 

of  the  level  of  significance  and  the  ratios  o^/o^,  Og/o^  of  the  variance 
components . 


-Ul- 


-.v-'  • 

.  Cj  1—  • ..  A 


' .  ^  ^ 


C.  ■<'.  w'-  Ji".  -  ■>'.  .V 1  -V. 


■  A-g-t. 


Power  values  associated  with  the  exact  test  statistic  F  were  computed  at 


the  a  =  . 10  level  of  significance.  The  upper  probability  values  of  the  F- 
distribution  (with  f  and  £2  degrees  of  freedom)  in  (3.5)  were  obtained  by 
using  the  IMSL  (International  Mathematical  and  Statistical  Libraries)  MDFDRE 
Subroutine  which  allows  fractional  degrees  of  freedom.  The  same  designs  and 
combinations  of  variance  components  as  the  ones  used  by  Tan  and  Cheng  (1984) 
were  considered  here.  The  results  are  given  in  Tables  2  and  3  for  designs  1 
and  2,  respectively.  From  these  tables  it  can  be  seen  tliat  the  exact  test  is 


more  efficient  than  the  approximate  tests  based  on  the  statistics  Fp  F2 ,  F^, 
and  F^.  Furthermore,  the  exact  test  has  the  advantage  I'.hat  its  critical  value 


does  not  depend  on  the  values  of  the  unknown  nuisance  parameters  a?  and  o^. 

Be 

This  property  is  not  shared  by  the  approximate  tests  which  may  require  good 


estimates  of  and  in  order  for  their  results  to  be  ri-llabie  (see  Tan  and 
p  e 


Cheng  1984,  p.  194). 


4 .  Concl I iding  remar  1^ 

The  vector  defined  in  (2.25)  is  the  key  to  the  ('onst  met  ion  of  the 
exact  test  proposed  in  Section  2.  It  is  a  linear  combination  of  the  random 
vectors  and  C';^,  where  ^  is  a  linear  transform  of  tlu-  vector  ^  of  response 
means,  and  C^^  makes  up  a  portion  of  the  residual  sum  of  sejuares  (see  formula 
2.16).  The  power  study  in  Section  3  clearly  indicates  that  the  exact  test  can 
be  at  least  as  efficient  as  the  other  approximate  tests. 

A  procedure  similar  to  the  one  desciibed  in  this  jiaper  was  effectively 
used  to  obtain  exact  tests  concerning  the  main  effects'  variance  components  in 
an  unbalanced  random  two-way  model  with  interai't  ion .  Details  of  that 
procedure  are  given  in  Khuri  (1985). 


Table  1 

Two  unbalanced  nested  designs 

for  the  model  in  (2.1) 

Design  1 

Design  2 

a=5 

a=4 

b^=2 

^1~'>  ^2=2 

b3=3,  b^=4 

nii=ni2=4 

nii=4 

"2r"3r"4r"5r^ 

"2r"22"^ 

"3l""32""33"^ 

"4r"42""43""44"^ 

Table  2 


Power  values  of  the  exact  and  approximate  tests  for  Design  1 
at  the  a  =  .10  level  of  significance. 


1 


.5 

0^  =  1.0 

[  T 

2  =  3. 

.0 

a 

a 

_  “  _ 

2 

^1 

•^2 

Exact 

^2 

^4 

Exact 

‘^l 

^2 

^4 

Exact 

1.0 

1. 

.0 

.099 

.097 

.098 

.118 

.099 

.097 

.094 

.134 

.  102 

.113 

.101 

.184 

l.O 

2, 

.0 

.099 

.099 

.097 

.114 

.098 

.101 

.095 

.126 

.098 

.129 

.097 

.167 

1.0 

3. 

.0 

.099 

.101 

.096 

.111 

.097 

.105 

.094 

.  121 

.097 

.136 

.096 

.156 

2.0 

1. 

,0 

.097 

w,.l53 

.099 

.111 

.097 

.154 

.099 

.121 

.098 

.164 

.10 

.156 

2.0 

2. 

.0 

.  10 

.098 

.099 

.109 

.099 

.097 

.098 

.118 

.10 

.098 

.097 

.148 

2.0 

3. 

.0 

.102 

.066 

.098 

.108 

.101 

.064 

.097 

.  1 16 

.101 

.061 

.095 

.142 

3.0 

1  , 

.0 

.  10 

.099 

.099 

.  108 

.10 

.099 

.099 

.116 

.  101 

.098 

.10 

.142 

3.0 

2, 

.0 

.10 

.099 

.099 

.  107 

.10 

.098 

.099 

.114 

.10 

.098 

.098 

.137 

3.0 

3, 

.0 

.  10 

.099 

.099 

.106 

.098 

.098 

.098 

.112 

.099 

.097 

.096 

.134 

Table  3 


Power  values  of  the  exact  and  approximate  tests  for  Design  2 
at  the  a  =  .10  level  of  significance 


_°L- 

1.0 

l  .0 

.150 

1.0 

2.0 

.139 

1.0 

3.0 

.130 

2.0 

1.0 

.127 

2.0 

2.0 

.123 

2.0 

3.0 

.120 

3.0 

1  .0 

.118 

3.0 

2.0 

.115 

3.0 

3.0 

.115 

159 

.141 

131 

,131 

126 

119 

118 


Exact 


147 

,135 

127 

125 

122 

118 

117 


.212 

.174 

.155 

.174 

.155 

.143 

.155 

.143 

.136 


°a  =  ^ 

.0 

=  3. 

0 

^'l 

^2 

^3 

Exact 

■=■1 

••’2 

^3 

Exact 

211 

.229 

.203 

.316 

.440 

.473 

.424 

.577 

185 

.194 

.177 

.248 

.374 

.427 

.353 

.472 

166 

.173 

.159 

.212 

.323 

.347 

.302 

.402 

158 

.166 

.155 

.248 

.301 

.323 

.294 

.472 

149 

.158 

.147 

.212 

.248 

.296 

.263 

.402 

143 

.145 

.139 

.189 

.251 

.267 

.238 

.353 

138 

.  140 

.137 

.212 

.236 

.245 

.231 

.402 

134 

-138 

.132 

.189 

.221 

.2  36 

.215 

.353 

131 

.136 

.130 

.174 

.210 

.227 

.203 

.316 

117  .115 
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